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An amplitude analysis of the 7r°7r° system produced in radiative J/ip decays is presented. In 
particular, a piecewise function that describes the dynamics of the 7r°7r° system is determined as a 
function of M n 0^0 from an analysis of the (1.311 ±0.011) x 10 9 J/ip decays collected by the BESIII 
detector. The goal of this analysis is to provide a description of the scalar and tensor components of 
the 7r°7r° system while making minimal assumptions about the properties or number of poles in the 
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amplitude. Such a model-independent description allows one to integrate these results with other 
related results from complementary reactions in the development of phenomenological models, which 
can then be used to directly fit experimental data to obtain parameters of interest. The branching 
fraction of J/ip —» 77r°7r° is determined to be (1.15±0.05) x 10 3 , where the uncertainty is systematic 
only and the statistical uncertainty is negligible. 

PACS numbers: 11.80.Et, 12.39.Mk, 13.20.Gd, 14.40.Be 


I. INTRODUCTION 

While the Standard Model of particle physics has 
yielded remarkable successes, the connection between 
the quantum chromodynamics (QCD) and the complex 
structure of hadron dynamics remains elusive. The light 
isoscalar scalar meson spectrum ( I G J PC = 0 + 0 ++ ), for 
example, remains relatively poorly understood despite 
many years of investigation. This lack of understand¬ 
ing is due in part to the presence of broad, overlapping 
states, which are poorly described by the most accessible 
analytical methods (see the “Note on scalar mesons below 
2 GeV” in the PDG) 0 . The PDG reports eight 0+0++ 
mesons, which have widths between 100 and 450 MeV. 
Several of these states, including the / o (1370), are char¬ 
acterized in the PDG only by ranges of values for their 
masses and widths. 

Knowledge of the low mass scalar meson spectrum is 
important for several reasons. In particular, the lightest 
glueball state is expected to have scalar quantum num¬ 
bers [1^5]. The existence of such a state is an excellent 
test of QCD. Experimental observation of a glueball state 
would provide evidence that gluon self-interactions can 
generate a massive meson. Unfortunately, glueballs may 
mix with conventional quark bound states, making the 
identification of glueball states experimentally challeng¬ 
ing. The low mass scalar meson spectrum is also of inter¬ 
est in probing the fundamental interactions of hadrons in 
that it allows for testing of Chiral Perturbation Theory 
to one loop [b]. 

The scalar meson spectrum has been studied in many 
reactions, including ttN scattering Q, pp annihilation df, 
central hadronic production Q, decays of the ip' jioj . 
J/ip fill-ill . B E, D [15| . and K [Til mesons, 77 for¬ 
mation [171 ] and (p radiative decays [181 ] . In particular, a 
coupled channel analysis using the K-matrix formalism 
has been performed using data from pion production, 
pp and np annihilation, and 7 T 7 t scattering [191 ]. Similar 
investigations would benefit from the inclusion of data 
from radiative J/ip decays, which provide a complemen¬ 
tary source of hadronic production. 

An attractive feature of a study of the two pseu¬ 
doscalar spectrum in radiative J/ip decays is the rela¬ 
tive simplicity of the amplitude analysis. Conservation of 
parity in strong and electromagnetic interactions, along 
with the conservation of angular momentum, restricts the 
quantum numbers of the pseudoscalar-pseudoscalar pair. 
Only amplitudes with even angular momentum and posi¬ 
tive parity and charge conjugation quantum numbers are 
accessible ( J PC = 0 ++ ,2 ++ ,4 ++ , etc). Initial studies 


suggest that only the 0 ++ and 2 ++ amplitudes are sig¬ 
nificant in radiative J/ip decays to 7r°7r°. The neutral 
channel (7r°7r°) is of particular interest due to the lack 
of sizable backgrounds like pn, which present a challenge 
for an analysis of the charged channel ( 7 r + 7 r“) ( 20 j |. 

Radiative J/ip decays to 7 r + 7 r~ have been analyzed pre¬ 
viously by the Marklll [2l| , DM2 [22j , and BES H3 ex¬ 
periments. Decays to 7r°7r° were also studied at Crystal 
Ball [24| and BES (25j, but these analyses were severely 
limited by statistics, particularly for the higher mass 
states. Each of these analyses reported evidence for 
the /2(1270) and some possible additional states near 
1.710 GeV/c 2 and 2.050 GeV/c 2 . More recently, the BE- 
SII experiment studied these channels and implemented 
a partial wave analysis [20| . Prominent features in the re¬ 
sults include the /2(1270), /o(1500), and /o(1710). How¬ 
ever, this analysis, like its predecessors, was limited by 
complications from large backgrounds and low statistics. 
Due to statistical limitations, the tt 0 tt 0 channel was used 
only as a cross check on the analysis of the charged chan¬ 
nel. 

Historically, amplitude analyses like that in Ref. poj 
have relied on modeling the s-dependence of the 7 T 7 t inter¬ 
action, where s is the invariant mass squared of the two 
pions, as a coherent sum of resonances, each described 
by a Breit-Wigner function. In doing so, a model is built 
whose parameters are resonance properties, e.g. masses, 
widths and branching fractions. A correspondence ex¬ 
ists between these properties and the residues and poles 
of the 7 T 7 r scattering amplitude in the complex s plane; 
however, this correspondence is only valid in the limit 
of an isolated narrow resonance that is far from open 
thresholds ( cf. Ref. 0 )- For regions containing mul¬ 
tiple overlapping resonances with large widths and the 
presence of thresholds, all of which occur in the 0 ++ 7 T 7 t 
spectrum, an amplitude constructed from a sum of Breit- 
Wigner functions becomes an approximation. While such 
an approximation provides a practical and controlled way 
to parameterize the data - additional resonances can be 
added to the sum until an adequate fit is achieved ■ it is 
unknown how well it maintains the correspondence be¬ 
tween Breit- Wigner parameters and the analytic struc¬ 
ture of the 7 T 7 r amplitude that one seeks to study, i.e., 
the fundamental strong interaction physics. Often sta¬ 
tistical precision, a lack of complementary constraining 
data, or a limited availability of models leaves the simple 
Briet-Wigner sum as a necessary but untested assump¬ 
tion in analyses, thereby rendering the numerical result 
only useful in the context of that assumption. In the 
context of this paper we refer to the Breit-Wigner sum 
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as a “mass dependent fit”, that is, the model used to fit 
the data has an assumed s dependence. 

In this analysis we exploit the statistical precision pro¬ 
vided by (1.311 ± 0.011) x 10 9 J/ip decays collected with 
the BESIII detector [26|, [13] to measure the components 
of the 7 T 7 r amplitude independently for many small re¬ 
gions of 7 T 7 T invariant mass, which allows one to con¬ 
struct a piecewise complex function from the measure¬ 
ments that describes the s- dependence of the 7 T 7 t dy¬ 
namics. Such a construction makes minimal assumptions 
about the s-dependence of the 7 T 7 t interaction. We refer 
to this approach in the context of the paper as a “mass 
independent fit”. 

The mass independent approach has some drawbacks. 
First, due to the large number of bins, one is left with 
a set of about a thousand parameters that describe the 
amplitudes with no single parameter tied to an individual 
resonance of interest. Second, mathematical ambiguities 
result in multiple sets of optimal parameters in each mass 
region. If only J = 0 and J = 2 resonances are signif¬ 
icant, there are two ambiguous solutions. However, in 
general, if one includes J > 4 the number of ambiguous 
solutions increases resulting in multiple allowed piecewise 
functions. Finally, in order to make the results practi¬ 
cally manageable for subsequent analysis, the assump¬ 
tion of Gaussian errors must be made - an assumption 
that cannot be validated in general. Similar limitations 
are present in other analyses of this type, e.g., Ref. Q. 
In spite of these limitations, which are discussed further 
in Appendices [B] and [C] the results of the mass indepen¬ 
dent amplitude analysis presented here represent a mea¬ 
surement of 7 T 7 T dynamics in radiative J/ij) decays that 
minimizes experimental artifacts and potential system¬ 
atic biases due to theoretical assumptions. The results 
are presented with the intent of motivating the devel¬ 
opment of dynamical models with reaction independent 
parameters that can subsequently be optimized using ex¬ 
perimental data. All pertinent information for the use of 
these results in the study of pseudoscalar-pseudoscalar 
dynamics is included in the supplemental materials (Ap¬ 
pendix 0 . 


II. THE BESIII DETECTOR 

The Beijing Spectrometer (BESIII) is a general- 
purpose, hermetic detector located at the Beijing 
Electron-Positron Collider (BEPCII) in Beijing, China. 
BESIII and BEPCII represent major upgrades to the BE- 
SII detector and BEPC accelerator. The physics goals 
of the BESIII experiment cover a broad research pro¬ 
gram including charmonium physics, charm physics, light 
hadron spectroscopy and r physics, as well as searches 
for physics beyond the standard model. The detector 
is described in detail elsewhere [28j . A brief description 
follows. 

The BESIII detector consists of five primary compo¬ 
nents working in conjunction to facilitate the reconstruc¬ 


tion of events. A superconducting solenoid magnet pro¬ 
vides a uniform magnetic field within the detector. The 
field strength was 1.0 T during data collection in 2009, 
but was reduced to 0.9 T during the 2012 running period. 
Charged particle tracking is performed with a helium-gas 
based multilayer drift chamber (MDC). The momentum 
resolution of the MDC is expected to be better than 0.5% 
at 1 GeV/c, while the expected dE/dx resolution is 6 %. 
With a timing resolution of 80 ps (110 ps) in the barrel 
(endcap), a plastic scintillator time-of-flight (TOF) de¬ 
tector is useful for particle identification. The energies of 
electromagnetic showers are determined using informa¬ 
tion from the electromagnetic calorimeter (EMC). The 
EMC consists of 6240 CsI(Tl) crystals arranged in one 
barrel and two endcap sections. With an angular cov¬ 
erage of about 93% of 47 t, the EMC provides an energy 
resolution of 2.5% (5%) at 1.0 GeV and a position resolu¬ 
tion of 6 mm (9 mm) in the barrel (endcap). Finally, par¬ 
ticles that escape these detectors travel through a muon 
chamber system (MUC), which provides additional infor¬ 
mation on the identity of particles. The MUC provides 
2 cm position resolution for muons and covers 89% of 
47 t. Muons with momenta over 0.5 GeV are detected 
with an efficiency greater than 90%. The efficiency of 
pions reaching the MUC is about 10% at this energy. 

Selection criteria and background estimations are stud¬ 
ied using a GEANT4 Monte Carlo (MC) simulation. The 
BESIII Object Oriented Simulation Tool (boost) [29| 
provides a description of the geometry, material compo¬ 
sition, and detector response of the BESIII detector. The 
MC generator KKMC [30j is used for the production of 
J/ip mesons by e + e~ annihilation, while BESEVTGEN |3l) 
is used to generate the known decays of the J/t/i accord¬ 
ing to the world average values from the PDG pj. The 
unknown portion of the J/ip decay spectrum is generated 
with the Lundcharm model [§3 |. 


III. EVENT SELECTION 

In order to be included in the amplitude analysis, an 
event must have at least five photon candidates and no 
charged track candidates. Any photon detected in the 
barrel (endcap) portion of the EMC must have an en¬ 
ergy of at least 25 (50) MeV. Four of the five photons are 
grouped into two pairs that may each originate from a 7r° 
decay. The invariant mass of any photon pair associated 
with a 7T° must fall within 13 MeV/c 2 of the 7r° mass. 
A 6 C kinematic fit is performed on each permutation of 
photons to the final state 77 r° 7 r°. This includes a con¬ 
straint on the four-momentum of the final state to that 
of the initial J/' i/j (4C) and an additional constraint (1C) 
on each photon pair to have an invariant mass equal to 
that of a 7T°. 

Significant backgrounds in this channel include J/if) de¬ 
cays to 777 (rj —> 7T 0 7r 0 7r°) and 777 ' ( r]' —> r)n°n 0 ; 77 — > 77 ). 
Restricting the y 2 from the 6 C kinematic fit is an ef¬ 
fective means of reducing the backgrounds of this type. 
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Events with a 7r°7r° invariant mass, M % » T o, below KK 
threshold (the region in which these backgrounds are sig¬ 
nificant) must have a y 2 less than 20. Events above KK 
threshold need only have a y 2 less than 60. To reduce 
the background from J/i/> decays to con 0 (w — t 77r°), the 
invariant mass of each jit 0 pair is required to be at least 
50 MeV/c 2 away from the u> mass [l|. Finally, in order 
to reduce the misreconstructed background arising from 
pairing the radiated photon with another photon in the 
event to form a 7r°, the invariant mass of the radiated 
photon paired with any n° daughter photon is required 
to be greater than 0.15 GeV/c 2 . 

If more than one permutation of five photons in an 
event satisfy these selection criteria, only the permuta¬ 
tion with the minimum y 2 from the 6C kinematic fit is 
retained. After all event selection criteria are applied, 
the number of events remaining in the data sample is 
442,562. MC studies indicate that the remaining back¬ 
grounds exist at a level of about 1.8% of the size of the 
total sample. Table Q] lists the major backgrounds. 

Backgrounds from J/iJ) decays to 7 r]C) are well un¬ 
derstood and are studied with an exclusive MC sample, 
which is generated according to the PDG branching frac¬ 
tions for these reactions. Other backgrounds are studied 
using an inclusive MC sample generated using BESEVT- 
GEN, with the exception of the misreconstructed back¬ 
ground, which is studied using an exclusive MC sample 
that resembles the data. The latter MC sample was gen¬ 
erated using a set of Breit-Wigner resonances with cou¬ 
plings determined from a mass dependent fit to the data 
sample. The °tv° spectrum after all selection criteria 
have been applied is shown in Fig.[Tj The reconstruction 
efficiency is determined to be 28.7%, according to the re¬ 
sults of the mass independent amplitude analysis. Con¬ 
tinuum backgrounds are investigated with a data sample 
collected at a center of mass energy of 3.080 GeV. The 
continuum backgrounds are scaled by luminosity and a 
correction factor for the difference in cross section as a 
function of center of mass energy. When scaled by lumi¬ 
nosity, only 3,632 events, which represents approximately 
0.8% of the signal, survive after all signal isolation re¬ 
quirements. 


IV. AMPLITUDE ANALYSIS 

A. General Formalism 

The results of the mass independent amplitude anal¬ 
ysis of the 7r 0 7r° system are obtained from a series of 
unbinned extended maximum likelihood fits. The ampli¬ 
tudes for radiative J/if) decays to 7r°7r° are constructed 
in the radiative multipole basis, as described in detail in 
Appendix [A] 

Let C/ M ’ A t represent the amplitude for radiative J/ip 
decays to 7r 0 7r°, 

(1) 


TABLE I. The number of events remaining after all selection 
criteria for each of a number of background reactions is shown 
in the right column. The backgrounds are broken into three 
groups. The first group contains the signal mimicking de¬ 
cays. The second lists the remaining backgrounds from J/i/j 
decays to 7 r/('), while the third group lists a few additional 
backgrounds. The backgrounds explicitly listed here represent 
about 93% of the total background according to the MC sam¬ 
ples. The misreconstructed background includes those events 
in which one of the daughter photons from a 7 r° decay is taken 
as the radiated photon. 


Decay channel 

Number of events 

J/t/j —> 77 r° 7 r° (data) 

442,562 

e + e —>■ 77 r° 7 r° (continuum) 

3,632 

J/ip —¥ bin 0 -, bi —> 77 T° 

1,606 

J/ip —> uin°-u> —» 77 T° 

865 

J/i/j —>■ pn°; p —> 77 T° 

778 

Misreconstructed background 

608 

J/ip — > 77 ; rj —> 37t° 

903 

J/ip —> 7 r/; r] —> rin°n°; r/ —» 77 

377 

J/ip —> a) 7 r° 7 r°;cj —» -yn 0 

775 

J/tf) —¥ &i7T°; bi —> ujn°-,uj —> 777 ° 

578 

J/ip —» w? 7 ;cj —» 77T° 

409 

J/ip — > 6l)/ 2 (1270); 00 — > 77 T° 

299 

J/i/j —> 7 7 c ;7c —► 'yn°n 0 orn 0 n°n 0 

255 

Other backgrounds 

507 

Total Background (MC) 

7,960 


where x = {0 7 , </> 7 , 8<p n } is the position in phase space, 
s = M^o^o is the invariant mass squared of the 7r°7r° pair, 
M is the polarization of the and A 7 is the helicity 
of the radiated photon. For the reaction under study the 
possible values of both M and A 7 are ±1. The amplitude 
may be factorized into a piece that contains the radiative 
transition of the J/ , 0 to an intermediate state X and a 
piece that contains the QCD dynamics 

U M ^{x,s)= £ (* 0 A H Q CD\X jt jJ 

( 2 ) 

where j is the angular momentum of the intermediate 
state and J 7 indexes the radiative multipole transitions. 
The sum over X includes any pseudoscalar-pseudoscalar 
final states (rnr, KK , etc) that may rescatter into 7r°7r°. 
We assume that the contribution of the 47 t final state to 
this sum is negligible, with the result that rescattering 
effects become important only above the KK threshold. 

The amplitude in Eq. © may be further factorized by 


U M ’ X ' 1 (x,s) = ( r yTr 0 TT°\H\J/ip) 
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FIG. 1. The M n o n o spectrum after all selection criteria have been applied. The black markers represent the data, while the 
histograms depict the backgrounds according to the MC samples. The signal (white) and misreconstructed background (pink) 
are determined from an exclusive MC sample that resembles the data. The other backgrounds are determined from an inclusive 
MC sample (see Tablc[I]). The components of the stacked histogram from bottom up are unspecified backgrounds, LL>7r°7r°, bin 0 , 
7 r/C), ujn °, the misreconstructed background, and the signal. 

pulling out the angular distributions, 

u M ^{x,s)= Y 

(3) 

where gj,j x{s) is the coupling for the radiative decay 
to intermediate state X. The functions ©j / ’ A ' T (0^, </>„•) 
and (0 7 , </> 7 ) contain the angular dependence of the 

decay of the X to 7r°7r° and the radiative J/V ) decay, re¬ 
spectively. The part of the amplitude that describes the 
7r°7r° dynamics is the complex function Tj t x{s), which is 
of greatest interest for this study. However, this func¬ 
tion cannot be separated from the coupling gj t j x(s)- 
Instead the product is measured according to 

v j,Jy( s ) ~y^gj,Jv,x(s)T jiX (s)- (4) 

This product will be called the coupling to the state with 
characteristics j, J 7 . Note here that, if rescattering ef¬ 
fects are assumed to be minimal (the only possible X is 
7T7r), all amplitudes with the same j have the same phase. 

The effect of rescattering is to break the factorizability 


of Eq. (|i)l . Finally, the amplitude may be written 

U M ^ (x, s) = Y Vj,J, (s)<c 7 (*), (5) 

where ( x ) contains the piece of the amplitude that 

describes the angular distributions and is determined by 
the kinematics of an event. 

Any amplitude with total angular momentum greater 
than zero will have three components (the 0 ++ amplitude 
has only an El component). Thus, three 2 ++ amplitudes, 
relating to El, M2, and E3 radiative transitions, are in¬ 
cluded in the analysis. While any amplitude with even 
total angular momentum and positive parity and charge 
conjugation is accessible for this decay, studies show that 
the 4 ++ amplitude is not significant in this region. In par¬ 
ticular, no set of four continuous 15 MeV/c 2 bins yield a 
difference in —2 In A greater than 28.8 units, which corre¬ 
sponds to a five sigma difference, under the inclusion of 
a 4 ++ amplitude. As no narrow spin-4 states are known, 
this suggests that only the 0 ++ and 2 ++ amplitudes are 
significant. The systematic uncertainty due to ignoring 
a 4 ++ amplitude that may exist in the data is described 
below in Sec. IV Cl 
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B. Parameterization 


C. Background subtraction 


The dynamical function in Eq. 0 may be parameter¬ 
ized in various ways. A common parameterization, dis¬ 
cussed in the introduction, is a sum of interfering Breit- 
Wigner functions, 

v j,J-y ( s ) = Yl k i’ J i’PBWj, Jy ,p(s), (6) 

0 


where BWjj p(s) represents a Breit-Wigner function 
with characteristics (mass and width) [3 and strength 

3- 

To avoid making such a strong model dependent as¬ 
sumption, we choose to bin the data sample as a function 
of M n o n o and to assume that the part of the amplitude 
that describes the dynamical function is constant over a 
small range of s, 

U M ^ (x, s) = J2 (a?). (7) 

3iJ 7 


For the scenario posed in Eq. ©, the couplings may 
be taken as the free parameters of an extended maximum 
likelihood fit in each bin of M n o T o. It is then possible to 
extract a table of complex numbers (the free parameters 
in each bin) that describe the dynamical function of the 
7 r° 7 T 0 interaction. 

The intensity function, I(x ), which represents the den¬ 
sity of events at some position in phase space x, is given 

by 


!(x)= E 


M.A-y 


3iJ'y 


(x) 


( 8 ) 


The mass independent amplitude analysis treats each 
event in the data sample as a signal event. For a clean 
sample, the effect of remaining backgrounds should be 
small relative to the statistical errors on the amplitudes. 
However, the backgrounds from J/xj) decays to 'yrjO in¬ 
troduce a challenge. Both of these backgrounds peak 
in the low mass region near interesting structures. The 
background from J /ip decays to 7 ?; lies in the region of 
the /o(500), which is of particular interest for its impor¬ 
tance to Chiral Perturbation Theory [l], [33|. The 777 ' 
background peaks near the /o(980), which is also of par¬ 
ticular interest due to its strong couplin g to I\K and its 
implications for a scalar meson nonet [341 ]. Therefore, 
the effect of these backgrounds is removed by using a 
background subtraction method. 

If a data sample is entirely free of backgrounds, the 
likelihood function is constructed as 

AT sig 

data 

m = n ( 10 ) 

i—1 

where /(£]£) is the probability density function (pdf) to 
observe an event with a particular set of kinematics x and 
parameters £ = {V^j }. The total number of parameters 
in the mass independent analysis is 1,178 (seven times 
the number of bins above KK threshold and five times 
the number of bins below KK threshold). The number 
of events in the pure data sample is given by N^ ta . 

Now, the likelihood may be written 

v slg v bkg JV bkg 

data data data 

^(D = n /(^i£) n /( f j 10 n ( n ) 

2 = 1 j =1 k—1 


The incoherent sum includes the observables of the re¬ 
action (which are not measured). For the reaction un¬ 
der study, the observables are the polarization of the 
J/il>, M = ±1, and the helicity of the radiated photon, 
Xy = ±1. The free parameters are constrained to be the 
same in each of the four pieces of the incoherent sum. 

In the figures and supplemental results that follow, the 
intensity of the amplitude in each bin is reported as a 
number of events corrected for acceptance and detector 
efficiency. That is, for the bin of M^o^o indexed by k and 
bounded by Sk and s^+i (the boundaries in s of the bin) 
we report, for each amplitude indexed by j and J 7 , the 
quantity 




dx. 


( 9 ) 


In practice, we absorb the size of phase space into the 
fit parameters. In doing so we fit for parameters V^j 
which are the V^j scaled by the square root of the size 
of phase space in bin k. 


where an additional likelihood, which describes the re¬ 
action for background events, has been multiplied and 
divided. Consider now a more realistic data sample that 
consists not only of signal events, but also contains some 
number of background events, Ay^,. Then the product 
of the first two factors of Eq. (fill) are simply the likeli¬ 
hood for the entire (contaminated) data sample, but the 
overall likelihood represents only that of the pure signal 
since the background likelihood has been divided. For a 
given data set, any backgrounds remaining after selection 
criteria have been applied are difficult to distinguish from 
the true signal. Rather than using the true background 
to determine the background likelihood, it is therefore 
necessary to approximate it with an exclusive MC sam¬ 
ple. That is, 

v bkg AT bkg 

-"data JV MC 

n fmr 1 ^ n /(*(12) 

2=1 2=1 

where the weight, Wi , is necessary for scaling purposes. 
For example, if the MC sample is twice the size of the 





expected background, a weight factor of 0.5 is necessary. 
Finally, the likelihood function may be written 


-fVdata "MC 

m = n /(*id n ( i3 ) 

*=1 3 = 1 

In practice, this likelihood distribution is multiplied by 
a Poisson distribution for the extended maximum likeli¬ 
hood fits such that 


m 


e ^ Ar ‘ lat 

-^data! 


a 


'data iV MC 

n /(* io n 


(14) 


An exclusive MC sample for the backgrounds due to 
J/if decays to 777 (') is generated according to the branch¬ 
ing fractions given by the PDG pj. This MC sample is 
required to pass all of the selection criteria that are ap¬ 
plied to the data sample. Any events that remain are 
included in the unbinned extended maximum likelihood 
fit with a negative weight (—Wj = — 1 in Eq. (fl 3 l) ~l. In this 
way, the inclusion of the MC sample in the fit approxi¬ 
mately cancels the effect of any remaining backgrounds 
of the same type in the data sample. 


D. Ambiguities 

Another challenge to the amplitude analysis is the 
presence of ambiguities. Since the intensity function, 
which is fit to the data, is constructed from a sum of 
absolute squares, it is possible to identify multiple sets of 
amplitudes which give identical values for the total inten¬ 
sity. In this way, multiple solutions may give comparable 
values of —21nL for a particular fit. For this particular 
analysis, two types of ambiguities are present. Trivial 
ambiguities arise due to the possibility of the overall am¬ 
plitude in each bin to be rotated by 7 r or to be reflected 
over the real axis in the complex plane. These may be 
partially addressed by applying a phase convention to the 
results of the fits. Non-trivial ambiguities arise from the 
freedom of amplitudes with the same quantum numbers 
to have different phases. The non-trivial ambiguities rep¬ 
resent a greater challenge to the analysis and cannot be 
eliminated without introducing model dependencies. 

While it is not possible in principle to measure the 
absolute phase of the amplitudes, it is possible to study 
the relative phases of individual amplitudes. Therefore in 
each of the fits, one of the amplitudes (the 2 ++ El ampli¬ 
tude) is constrained to be real. The phase difference be¬ 
tween the other amplitudes and that which is constrained 
can then be determined in each mass bin. 

As mentioned above, a set of trivial ambiguities arises 
due to the possibility of the overall amplitude in each bin 
to be rotated by 7 r or to be reflected over the real axis 
in the complex plane. Each of these processes leave the 
intensity distribution unchanged. This issue is partially 
resolved by establishing a phase convention in which the 


amplitude that is constrained to be real is also con¬ 
strained to be positive. The remaining ambiguity is re¬ 
lated to the inability to determine the absolute phase. 
The phase of the total amplitude may change sign with¬ 
out inducing a change in the total intensity. Therefore, 
when a phase difference approaches zero, it is not pos¬ 
sible to determine if the phase difference should change 
sign. The amplitude analysis results are presented here 
with the arbitrary convention that the phase difference 
between the 0 ++ amplitude and the 2 ++ El amplitude is 
required to be positive. One may invert the sign of this 
phase difference in a given bin, but then all other phase 
differences in that bin must also be inverted. 

The presence of non-trivial ambiguities is attributed to 
rescattering effects, which allow for amplitudes with the 
same quantum numbers, J PC , to have different phases. 
The couplings, x(s)> hr Eq- © are real functions 
of s. Since the dynamical amplitude, Tj x{s), does not 
depend on J 7 , its phase is the same for each of the am¬ 
plitudes with the same J PC (in particular, the 2 ++ El, 
M2 and E3 amplitudes). However, if more than one in¬ 
termediate state, A', is present, differences between the 
couplings to these amplitudes may result in a phase dif¬ 
ference. Therefore, in the region above the KK threshold 
the 2 ++ amplitudes may have different phases. However, 
below KK threshold the phases of these amplitudes are 
constrained to be the same. That is, rescattering through 
47 t is assumed to be negligible. 

By writing out the angular dependence of the intensity 
function, it is possible to show that the freedom to have 
phase differences between the components of a given am¬ 
plitude (2 ++ El, M2, and E3, for example) generates an 
ambiguity in the intensity distribution. For this chan¬ 
nel and considering only 0 ++ and 2 ++ amplitudes, two 
non-trivial ambiguous solutions may be present in each 
bin above KK threshold. The knowledge of one solu¬ 
tion can be used to mathematically predict its ambiguous 
partner. In fact, some bins do not exhibit multiple so¬ 
lutions, but have a degenerate ambiguous pair. A study 
of these ambiguities (Appendix [Bj shows consistency be¬ 
tween the mathematically predicted and experimentally 
determined ambiguities. Both ambiguous solutions are 
presented, because it is impossible to know which rep¬ 
resent the physical solutions without making some addi¬ 
tional model dependent assumptions. If more than two 
solutions are found in a given bin, all solutions within 1 
unit of log likelihood from the best solution are compared 
to the predicted value derived from the best solution and 
only that which matches the prediction is accepted as the 
ambiguous partner. 


E. Results 

1. Amplitude intensities and phases 

The intensity for each amplitude as a function of M r o w o 
is plotted in Fig.[2j Each of the phase differences with re- 
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spect to the reference amplitude (2 ++ El), which is con¬ 
strained to be real, is plotted in Fig. [3] Above the KK 
threshold, two distinct sets of solutions are apparent in 
most bins as expected. The bins below about 0.6 GeV/c 2 
also contain multiple solutions, but with different likeli¬ 
hoods and are attributed to local minima in the likeli¬ 
hood function. The nominal solutions below 0.6 GeV/c 2 
are determined by requiring continuity in each intensity 
and phase difference as a function of M n o n o. Only sta¬ 
tistical errors are presented in the figures. 

It is apparent that the ambiguous sets of solutions in 
the nominal results are distinct in some regions, while 
they approach and possibly cross at other points. The 
most powerful discriminator of this effect is the phase 
difference between the El and M2 components of the 
2 ++ amplitude (see the middle plot of Fig. Re¬ 
gions in which the solutions may cross are apparent at 
0.99 GeV/c 2 , near 1.3 GeV/c 2 , and above 2.3 GeV/c 2 . 
Since the results in each bin are independent of their 
neighbor, it is not possible to identify two distinct, 
smooth solutions at these crossings. 


2. Discussion 

The results of the mass independent analysis exhibit 
significant structures in the 0 ++ amplitude just below 
1.5 GeV/c 2 and near 1.7 GeV/c 2 . This region is where 
one might expect to observe the the states /o(1370), 
/o(1500), and /o(1710) which are often cited as being 
mixtures of two scalar light quark states and a scalar 
glueball §15], [Si]. A definitive statement on the number 
and properties of the scattering amplitude poles in this 
region of the spectrum requires model-dependent fits to 
the data. The effectiveness of any such mo del-dependent 
study could be greatly enhanced by including similar 
data from the decay J/ij) —> 'yKK in an attempt to iso¬ 
late production features from partial widths to KK and 
7T7T final states. 

Additional structures are present in the 0 ++ amplitude 
below 0.6 GeV/c 2 and near 2.0 GeV/c 2 . It seems reason¬ 
able to interpret the former as the a (/o(500)). The latter 
could be attributed to the /o(2020). The presence of the 
four states below 2.1 GeV/c 2 would be consistent with 
the previous study of radiative J/ip decays to 7T7 t by BE- 
SII [2y] . Finally, the results presented here also suggest 
two possible additional structures in the 0 ++ spectrum 
that were not observed in Ref. [lo]. These include a struc¬ 
ture just below 1 GeV/c 2 , which may indicate an /o(980), 
but the enhancement in this region is quite small. There 
also appears to be some structure in the 0 ++ spectrum 
around 2.4 GeV/c 2 . 

In the 2 ++ amplitude, the results of this analysis in¬ 
dicate a dominant contribution from what appears to 
be the /2(1270), consistent with previous results [20j |. 
However, the remaining structure in the 2 ++ amplitude 
appears significantly different than that assumed in the 
model used to obtain the BESII results 0 . In particu¬ 


lar, the region between 1.5 and 2.0 GeV/c 2 was described 
in the BESII analysis with a relatively narrow / 2 (1810). 
One permutation of the nominal results (the red markers 
in Fig. H indicates that the structures in this region are 
much broader, while the other permutation (the black 
markers in Fig. suggests that there is very little con¬ 
tribution from any 2 ++ states in this region. 

The tensor spectrum near 2 GeV/c 2 is of interest in 
the search for a tensor glueball. Previous investiga¬ 
tions of the J/if) —> 77r°7T 0 channel reported evidence 
for a narrow (T ss 20 MeV) tensor glueball candidate, 
fj( 2230) [25(. While a mo del-dependent fit is required 
to place a limit on the production of such a state using 
these data, we note that based on the reported value of 
—f 7/j( 2230)) [23j, one would naively expect to 
observe a peak for the fj( 2230) with an integral that is 
of order 4 x 10 5 but concentrated only in roughly two 
bins of M(7T 0 7r°), corresponding to the full width of the 
fj( 2230). Such a structure seems difficult to accommo¬ 
date in the extracted 2 ++ amplitude. 

F. Branching fraction 

The results of the mass independent amplitude analy¬ 
sis allow for a measurement of the branching fraction of 
radiative J/i^ decays to 7r°7r°, which is determined ac¬ 
cording to: 

77r u 7r u ) = —l ---(15) 

e 7 N J/ii> 

where N^ n o n o is the number of acceptance corrected 
events, A^kg is the number of remaining background 
events, e 7 is an efficiency correction necessary to extrap¬ 
olate the 7r°7r° spectrum down to a radiative photon en¬ 
ergy of zero, and Nj/^ is the number of J/ip decays in 
the data. The number of acceptance corrected events 
is determined from the amplitude analysis by summing 
the total intensity from each M n o n o bin. The number 
of remaining background events is determined according 
to the inclusive and exclusive MC samples. The frac¬ 
tional background contamination in each bin i. l?bkg,i, 
is determined before acceptance correction. The number 
of background events is then determined by assuming 
i?bkg,i is constant after acceptance correction such that 
the number of background events in bin i, iVbkg.i, is given 
by the product of i?bkg,i and the number of acceptance 
corrected events in the same bin, N~^o n o Note that the 
backgrounds from to J/ij) decays to are removed 

during the fitting process and are not included in this 
factor. The efficiency correction factor, e 7 , is determined 
by calculating the fraction of phase space that is removed 
by applying the selection requirements on the energy of 
the radiative photon. This extrapolation increases the 
total number of events by 0.07%. Therefore, e 7 is taken 
to be 0.9993. 

The backgrounds remaining after event selection fall 
into three categories. The misreconstructed backgrounds 
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FIG. 2. The intensities for the (a) 0 ++ , (b) 2 ++ El, (c) 2 ++ M2 and (d) 2 ++ E3 amplitudes as a function of M n o n o for the 
nominal results. The solid black markers show the intensity calculated from one set of solutions, while the open red markers 
represent its ambiguous partner. Note that the intensity of the 2 ++ E3 amplitude is redundant for the two ambiguous solutions 
(see Appendix 11). Only statistical errors are presented. 


are determined from an exclusive MC sample that re¬ 
sembles the data. Events that remain in a continuum 
data sample taken at 3.080 GeV after selection criteria 
have been applied are also taken as a background. Fi¬ 
nally, the other remaining backgrounds are determined 


using the inclusive MC sample. Each of these back¬ 
grounds is scaled appropriately. In total, the acceptance 
corrected number of background events, -/Vbkg) is deter¬ 
mined to be 35,951. The number of radiative J/ip decays 
to 7r°7r°, N 11 ro^o, is determined to be 1,543,050 events. 
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FIG. 3. The phase differences relative to the reference amplitude (2 ++ El) for the (a) 0 ++ , (b) 2 ++ M2, and (c) 2 ++ E3 
amplitudes as a function of M n o^o for the nominal results. The solid black markers show the phase differences calculated from 
one set of solutions, while the open red markers represent the ambiguous partner solutions. An arbitrary phase convention is 
applied here in which the phase difference between the 0 ++ and 2 ++ El amplitudes is required to be positive. Only statistical 
errors are presented. 


The branching fraction for this decay is then determined 
to be (1.151 ± 0.002) x 10 -3 , where the error is statistical 
only. 


V. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties for the mass independent 
analysis include two types. First, the uncertainty due 
to the effect of backgrounds from J/ip decays to 7 ry(') 
are addressed by repeating the analysis and treating the 
background in a different manner. The second type of 
systematic uncertainty is that due to the overall normal¬ 
ization of the results. Sources of systematic uncertainties 
of this type include the photon detection efficiency, the 
total number of J/ip decays, the effect of various back¬ 
grounds, differences in the effect of the kinematic fit be¬ 
tween the data and MC samples and the effect of model 


dependencies. The uncertainty on the branching fraction 
of 7T° to 77 according to the PDG is 0.03% [l|, which 
is negligible in relation to the other sources of uncer¬ 
tainty. The systematic uncertainties are described below 
and summarized in Table ITT1 These uncertainties also ap¬ 
ply to the branching fraction measurement. Finally, sev¬ 
eral cross checks are also performed. 


A. J / tp —> 777 and J/ip —> 7 r\ Background 

Uncertainty 

The amplitude analysis is performed with the assump¬ 
tion that all backgrounds have been eliminated. Stud¬ 
ies using Monte Carlo simulation indicate this is a valid 
assumption for most of the M r o T o spectrum. However, 
significant backgrounds from J/ip decays to 777 and 777 ' 
exist in many mass bins below about 1 GeV/c 2 . Rather 
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than inflating the errors of these bins according to the un¬ 
certainty introduced by these backgrounds, which would 
not take into account the bin-to-bin correlations, a set 
of alternate results is presented in which the 'yrjC) back¬ 
grounds are not subtracted. 

The fraction of events in J/if decays to 7 rj(') that sur¬ 
vive the event selection criteria for the 77 r° 7 r° final state 
is very small (about 0.02%). Minor changes to the mod¬ 
eling of these decays may therefore have a large effect 
on the backgrounds. The difference between the nominal 
results and the alternate results, which treat the back¬ 
grounds differently, can be viewed as an estimator of the 
systematic error in the results due to these backgrounds. 

The distinctive feature of the alternate results is an 
enhancement in the 0 ++ intensity in the region below 
about 0.6 GeV/c 2 and near the 77 ' peak. This may be 
interpreted as the contribution of the events from J/•if 
decays to y^), which are being treated as signal events. 
A comparison of the 0“ l_+ amplitude for nominal results 
and the alternate results is presented in Fig. [I] The re¬ 
sults for the other amplitudes are consistent between the 
two methods. Any conclusion drawn from these data 
that is sensitive to choosing specifically the alternate or 
nominal results is not a robust conclusion. 


B. Uncertainties in the overall normalization 

1. Photon Detection Efficiency 

The primary source of systematic uncertainty for this 
analysis comes from the reconstruction of photons. To 
account for this uncertainty, the photon detection effi¬ 
ciency of the BESIII detector is studied using the so 
called tag and probe method on a sample of J/if de¬ 
cays to 7r + 7r _ 7T°, where the 7r° decays into two photons. 
One of these final state photons is reconstructed, along 
with the two charged tracks, while the other photon is 
left as a missing particle in the event. This information 
can then be used to determine the region in the detec¬ 
tor where the missing photon is expected. The photon 
detection efficiency is calculated by taking the ratio of 
the number of missing photons that are detected in this 
region to the number that are expected. The numbers of 
detected and expected photons are determined with fits 
to the two photon invariant mass distributions. 

The systematic uncertainty due to photon reconstruc¬ 
tion is determined by investigating the differences be¬ 
tween the photon detection efficiencies of the inclusive 
MC sample and that of the data sample. This difference 
is measured to be less than 1 . 0 %, which is taken to be 
the systematic uncertainty per photon. For the five pho¬ 
ton final state the overall uncertainty due to this effect 
is therefore taken to be 5.0%. 

An additional source of uncertainty, which is due to 
mismodelling of the photon detection efficiency as a func¬ 
tion of the angular and energy dependence of the radia¬ 
tive photon, was studied using the same channel. The 


phase space MC samples used for normalization in each 
bin of the mass independent amplitude analysis were 
modified to account for differences in the photon detec¬ 
tion efficiency between the data and inclusive MC sam¬ 
ples. The mass independent analysis was then repeated 
using the modified phase space MC samples. Neither the 
differences in angular nor energy dependence had a sig¬ 
nificant effect on the results of the analysis. The effects 
of mismodelling of this type are therefore taken to be 
negligible. 


2. Number of J/if 


The number oi J /if decays is determined from an anal¬ 
ysis of inclusive hadronic events 


Nj/i/j 


N se \ — N b g 


^trig 


A(2S) , 

Mata A 


(16) 


where N se \ represents the number of inclusive events re¬ 
maining after selection criteria have been applied and 
N bg is the number of background events estimated with a 
data sample collected at 3.080 GeV. The efficiency for the 
trigger is given by etrig, while e^ita^ the detection ef¬ 
ficiency for J/if inclusive decays determined from if(2S) 
decays to tt + tt~ J/if . Finally, / cor represents a correction 
factor to translate e<fata ^ to the efficiency for inclusive de¬ 
cays in which the J/if is produced at rest. To obtain N se \, 
at least two charged tracks are required for each event. 
Additionally, the momenta of these tracks and the visible 
energy of each event are restricted in order to eliminate 
Bhabha and di-muon events as well as beam gas inter¬ 
actions and virtual photon-photon collisions. The total 
number of J/if decays in the data sample according to 
Eq. (fT 6 l) is determined to be (1.311 ± 0.011) x 10 9 events, 
which results in an uncertainty of 0.8% [26l |27I| . 


3. Background Size 

According to the inclusive MC sample, the total num¬ 
ber of background events that contaminate the signal is 
about 1.5%. These do not include the misreconstructed 
backgrounds nor the backgrounds from J/if decays to 
7 ??('), both of which are addressed in a separate system¬ 
atic uncertainty. Additionally, backgrounds from non- 
J/if decays yield a contamination of approximately 0.8%. 
Conservative systematic uncertainties equal to 100% of 
the background contamination are attributed to each of 
the inclusive MC and continuum background types. 


4- Uncertainty in the acceptance corrected signal yield 

One of the largest remaining backgrounds after signal 
isolation and background subtraction is the signal mim¬ 
icking decay of J/if to W7r°, where the w decays to 77r°. 
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FIG. 4. A comparison of the (a) 0 ++ intensity and (b) phase difference relative to the 2 ++ El amplitude for the nominal results 
and the alternate results, in which the 7 r](') backgrounds have not been subtracted from the data. The solid black markers 
show the nominal results, while the red markers represent the alternate results. Only statistical errors are presented. 


The nominal method to address this background is to re¬ 
strict the 77 T° invariant mass to exclude the region within 
50 MeV/c 2 of the w mass. An alternative method is to in¬ 
clude an amplitude for the W 7 r° final state in the analysis. 
The results of this alternative method are quantitatively 
no different than the nominal results, suggesting that the 
exclusion method is an effective means of addressing the 
background from J/ip decays to uin 0 . The difference in 
the branching fraction using the signal yield for the alter¬ 
native method compared to the nominal method is about 
0 . 8 %. 

As discussed above, backgrounds due to J/^> decays 
to 7 riC) are addressed in the fitting procedure itself by 
adding an exclusive MC sample to the data, but with a 
negative weight. The systematic uncertainty do to this 
background is determined by using the data alone. In this 
way, contributions from these backgrounds are treated as 
signal and inflate the signal yield and background size 


in Eq. m- The difference in the branching fraction is 
0.03%, which is considered a negligible contribution to 
the systematic uncertainty. 

Differences in the effect of the 6 C kinematic fit on the 
data and MC samples may cause a systematic difference 
in the acceptance corrected signal yield. This effect was 
investigated by loosening the restriction on the y 2 from 
the 6 C kinematic fit. For events with a M n o n o above 
KK threshold, this restriction was relaxed from less than 
60 to be less than 125. Events with an invariant mass 
below KK threshold are required to have a y 2 less than 60 
rather than less than 20. The difference in the branching 
fraction for the results with the loosened y 2 cut relative 
to that of the nominal results is about 0 . 1 %. 

Another source of systematic uncertainty in the 
branching fraction is the difference between the nomi¬ 
nal results and those obtained by applying a model that 
describes the 7 T 7 t dynamics. To test this effect, a mass 
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dependent fit using interfering Breit-Wigner line shapes 
was performed. The difference in the branching fraction 
using the acceptance corrected yield of the mass depen¬ 
dent analysis compared to the nominal results is about 
0.3%. 

The effect of the remaining misreconstructed back¬ 
grounds on the results is studied by performing a closure 
test, in which the mass independent amplitude analysis 
is performed on an exclusive MC sample. This MC sam¬ 
ple was generated according to the results of a mass de¬ 
pendent amplitude analysis of the data and includes the 
proper angular distributions. After applying the same 
selection criteria that are applied to the data, the MC 
sample is passed through the mass independent analy¬ 
sis. This process is repeated after removing the remain¬ 
ing misreconstructed backgrounds from the sample. The 
difference in the branching fraction between these two 
methods is 0.01%. The effect of these backgrounds is 
therefore taken to be negligible. 

TABLE II. This table summarizes the systematic uncertain¬ 
ties (in %) for the branching fraction of radiative J/ip decays 

to 7r°7r°. 


Source 

J/tp —¥ 77r°7r° (%) 

Photon detection efficiency 

5.0 

Number of J /■ip 

0.8 

Inclusive MC backgrounds 

1.5 

Non- J/ip backgrounds 

0.8 

ojn° background 

0.8 

Kinematic fit X6C 

0.1 

Model dependent comparison 

0.3 

Total 

5.4 


C. 4 ++ amplitude 

As discussed above, the only 7r°7r° amplitudes that 
are accessible in radiative J/ip decays have even angu¬ 
lar momentum and positive parity and charge conjuga¬ 
tion quantum numbers. The mass independent analysis 
was performed under the assumption that only the 0 ++ 
and 2 ++ amplitudes are significant. To test this assump¬ 
tion, the analysis was repeated with the addition of a 
4 ++ amplitude. No significant contribution from a 4 ++ 
amplitude is apparent. 

To test the effect of a 4 ++ amplitude that may exist in 
the data and is ignored in the fit, an exclusive MC sample 
was generated using a model constructed from a sum of 
resonances each parameterized by a Breit-Wigner func¬ 
tion in a way that optimally reproduces the data. One 
of the resonances was an /4(2050), which was generated 
in each component of the 4 ++ amplitude. The relative 
size of the 4 ++ amplitude was determined from a mass 
dependent fit to the data, in which the 4 ++ amplitude 
contributed 0.43% to the overall intensity. A mass inde¬ 


pendent amplitude analysis, which did not include a 4 ++ 
amplitude, was then performed on this sample. The re¬ 
sults indicate that the intensities and phases for the 0 ++ 
and 2 ++ amplitudes deviate from the input parameters 
at the order of the statistical errors from the data sample 
in the region between 1.5 and 3.0 GeV/c 2 . Therefore, the 
systematic error due to the effect of ignoring a possible 
4 ++ amplitude is estimated to be of the same order as 
the statistical errors in the region from 1.5 to 3.0 GeV/c 2 . 

VI. CONCLUSIONS 

A mass independent amplitude analysis of the 7r°7r° 
system in radiative J/ip decays is presented. This anal¬ 
ysis uses the world’s largest data sample of its type, col¬ 
lected with the BESIII detector, to extract a piecewise 
function that describes the scalar and tensor 7T7 t ampli¬ 
tudes in this decay. While the analysis strategy employed 
to obtain results has complications, namely ambiguous 
solutions, a large number of parameters, and potential 
bias in subsequent analyses from non-Gaussian effects 
(see Appendix [Cj, it minimizes systematic bias arising 
from assumptions about 7T7 t dynamics, and, consequently, 
permits the development of dynamical models or param- 
eterizations for the data. 

In order to facilitate the development of models, the 
results of the mass independent analysis are presented in 
two ways. The intensities and phase differences for the 
amplitudes in the fit are presented here as a function of 
M n o^o. Additionally, the intensities and phases for each 
bin of o w o are given in supplemental materials (see 
Appendix [UJ. These results may be combined with those 
of similar reactions for a more comprehensive study of 
the light scalar meson spectrum. Finally, the branching 
fraction of radiative J/ip decays to 7r°7r° is measured to 
be (1.15 ± 0.05) x 10~ 3 , where the error is systematic 
only and the statistical error is negligible. This is the 
first measurement of this branching fraction. 
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Appendix A: Amplitudes 

The amplitude for radiative J/f/' decays to 7r°7r° can 
be determined in different bases depending on the infor¬ 
mation of interest. For example, in the lrelicity basis, 
the amplitude depends on the angular momentum and 
helicity of the 7r 0 7r° resonance as well as the angular mo¬ 
mentum and polarization of the J/ip. It is also possible to 
relate the amplitudes to radiative multipole transitions. 
Such a basis is useful because it may allow implementa¬ 
tion or testing of dynamical assumptions. For example, 
a model may suggest that the El radiative transition 
should dominate over the M2 transition. 

In the radiative multipole basis, the amplitude for ra¬ 
diative J/ip decays to 7r°7r° is given by 


U M ^{x,s) = Y, Nj,N 3 D j m ^ (tt + 0 7 , tt - fl 7 , 0)D^ o (^, 0„, 0)l 1 + ( lV 
j s J'Y > M 

(J 7 — A 7 ; jfi\ — A 7 ) -y^[c>A 7 ,i ]Vj ,J 7 ( S ) 


(Al) 


where the parity, total angular momentum, and helicity 
of the pair of pseudoscalars are given by P, j, and /z, 
respectively. The D functions are the familiar Wigner 
D-matrix elements. The angular momentum of the pho¬ 
ton, J 7 , is related to the nuclear radiative (El, M2, E3, 
etc.) transitions. Each amplitude is characterized by the 
angular momentum of the photon and the angular mo¬ 
mentum of the pseudoscalar pair. The possible values of 
J 7 are limited by the conservation of angular momentum. 
The helicity of the radiative photon is given by A 7 . The 
total angular momentum and polarization of the J/ij) are 

given by J and M, respectively. Finally, Nj = yj is 
a normalization factor. 

The angles (</> 7 , 0 7 ) are the azimuthal and polar angles 
of the photon in the rest frame of the J/ip, where the 
direction of the J/ip momentum defines the x-axis. The 
angles (<p v , 9 n ) are the azimuthal and polar angles of one 
7 r° in the rest frame of the 7r°7r° pair, with the -z axis 
along the direction of the photon momentum and the x- 
axis is defined by the direction perpendicular to the plane 
shared by the beam and the z-axis. 

Parity is a conserved quantity for strong and elec¬ 
tromagnetic interactions. Hence, for J/ip radiative de¬ 
cays, P = (—1) J must be positive. This means that the 
only intermediate states available have j p = 0 + ,2 + ,4 + , 
etc. Additionally, isospin conservation in strong inter¬ 
actions requires I G for the intermediate state to be 0 + 
(isoscalar). The complex function U,,j 7 (s) describes the 


7r°7r° production and decay dynamics. In order to min¬ 
imize the model dependence of the mass independent 
analysis, the dynamical amplitude is replaced by a (com¬ 
plex) free parameter in the unbinned extended maximum 
likelihood fit. Thus, the amplitude, in a region around s 
is given by 


U M ^ (x, s) = Y (f), (A2) 

31 J'Y 


where 


= N J-r N 3 D M,p- A 7 (P + hi'* - 0 7>°) 
11 + f—IV 

K, 0 (6rA,0)g l 2 ^ 

(</ 7 A 7 , jfi\J /-i A 7 ) 


(A3) 


V2 




and {j, J 7 } represents the unique amplitudes accessible 
for the given set of observables, {M, A 7 }. 


Appendix B: Ambiguities 

One of the challenges of amplitude analysis is the issue 
of ambiguous solutions, two solutions that give the same 
distribution (eg. Ref. Q). In this section, the ambiguous 
solutions for radiative J/'i/' decays to 7r°7r° are studied. 
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To determine the angular dependence of the ampli¬ 
tudes, it is necessary to write the decay amplitude 
( x ), which is given in Eq. JAD, explicitly as a func¬ 
tion of the angles (<^ 7 , 0 7 ) and 9 n ). The Clebsch 
Gordan factors in the amplitude restrict the signs of /x to 


be the same as that of A 7 . Thus, for j = 2 and A 7 = 1, 
only the values /i = 0 , 1 , 2 give non-zero amplitude con¬ 
tributions. It is also important to note that the Clebsch 
Gordan coefficients will change sign under A 7 —> — A 7 , 
but only for J 7 = 2. This will cancel the delta functions 
in the decay amplitude with the result 


= y : yyNj^N je - iM ^ + ^dy_ x ^ - 9,) x e-^4, 0 (<4)^#A 71 i + ^(-l) 7 ^ 1 ] (Bl) 




H,0\ 


V2 l 


r 


where the constants cf),’ 7 contain the Clebsch-Gordan 

Jlh 1 

coefficients. 

Recall that, for the Wigner small d-matrix elements, 
d i,±i( n — 9) = d\^(9) and d\ 0 (-K — = d i,o( d )- Then, 

d M^-\y ~ 6 ) = d M, A 7 - /i ( 6 ')- Also > note that the re¬ 
strictions on /x mean that the quantity /j, — A 7 = ±1,0. 
It is also useful to note that /x — A 7 = A 7 ,0, — A 7 , for 
/x = ±2, ±1,0 respectively. The usefulness of these fea¬ 
tures appears when one writes out the intensity for a 
given choice of M and A 7 . It is also useful to plug in 
the values for the constants, which are given in Table ITITI 
The intensity in bin a for a given choice of observables is 
then given by 

i a (x) = y 


M,\~y 


+ h2(9n)d 1 M x (6~/)e 


— iXry&TT |2 


(B2) 


where terms with the same angular dependencies have 
been grouped according to 

ho(Qn) = >/3Vo, 1 ± \j 7 ^(V2,l ± a/ 5V2,2 ± 2I / 2 I 3)do j o(^7r) 
hi(9„) = -^=( 3 ^ 2,1 ± V5V 2 , 2 - 4 V 2 , 3 )dl 0 (9„) 
h 2 (0 n ) = (3 V 2 ,i - V5V 2 ,2 ± V 2 , s )d\ o (0„) 

(B3) 

and the subscripts on the production amplitudes repre¬ 
sent the possible combinations of j and J 1 . The following 
calculations apply for each bin individually. 

TABLE III. The constant factors in Eq. (ED are given here. 




= 1 


AA 1 _ 

l 2,0 — 

,±H - 
l 2,1 — 


r 2 ,±i - ±, 

l 2,0 — ^ 


.2>±1 - 


2,±1 


u 2,2 


f3.il 

l 2,0 


1 n ^2 ^ 1/ 3K 

c 3.±l - 
l 2,2 ‘ 


The amplitudes for which M and A 7 have the same 
(opposite) sign, M = A 7 = ±1 (M = —A 7 = ±1) are 


related to each other by a sign change in the exponential 
factor. Note that the terms with a factor of d x M 0 will 
change sign under M —t —M and terms with a factor 
of 4 0 will change sign under A 7 —> — A 7 . Then, the 
intensity becomes 

I(x)= y |fto(^)di,i(0 7 )e ±i ^+/n(^)di,o(^) 

M=A 7 =il 

+ h2(e n )d 1 lt _ 1 (e i )e^\ 2 

+ y \ho(6„)d\ t _ 1 (0 1 )e ±i +* -/ii(^) 4 ,o(^) 


M=-A 7 =il 


+ h2(e n )dlMe^\ 2 . 


(B4) 


Note that the term with hi(9^) has changed sign in 
the opposite combination. The properties of small d 
functions, d J m , m (9) = (-1 ) m ~ m 'd J m m ,(9) = d J _ m _ m ,(9), 
have been used to write the incoherent pieces of the in¬ 
tensity in the same way. 

It is instructive to write the intensity function as 


I(x) =/o ± fi cos 2 9 1 ± -f 2 COS 2 


± i / 3 sin 20 7 cos <j)„ — i /4 cos 2 9 1 cos 2 


(B5) 


where 


fo — ± (^ 2 ) 2 ] ± (h\) 

fi = yh 0 ) 2 + (h 2 ) 2 ] - (hi) 2 

h = fi = (hoh 2 + h* 0 h 2 ) 

h = V2(-hoK - h*hi + h 2 h\ ± h* 2 hi). 


(B 6 ) 


Now, if a set of amplitude couplings, V, have been 
determined by fitting the intensity function in Eq. (IB5I) to 
the data, ambiguities would arise if an alternative set of 
couplings, V', would give the same angular dependence as 
the original set. In other words, the new set of amplitudes 
must give the same values for the /,; functions (/' = fi). 

Consider f 2 , which can be written as a linear combi¬ 
nation of two quadratic forms 


fi — y(\ho ± h 2 \ 2 — | ho — h 2 \ 2 ). 


(B7) 
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These quadratic forms are given by 


| h 0 ± h 2 1 2 = [cos 2 9 n (3ai q= a 3 ) + {b - a 1 ± a 3 )] 
x [cos 2 6^(3 a\ =F a*) + (b* - a\ ± a^)], 


(B8) 


where for simplicity the production coefficients have been 
combined into new variables given by 


b = v / 3Vo,i 
y/6, 


ai — —— (V^i + -\/5V2,2 + 2V2 t3 ) 


VS, 


a 2 = ~(SV 2 ,i + y/5V 2 , 2 - 4U 2 , 3 ) 


(B9) 


V6, 


a 3 = ^-{SV 2 ,i - VbV 2 , 2 + V 2 , 3 ). 


Since only the absolute square of each combination of 
ho and h 2 appears in the intensity, nontrivial ambiguous 
solutions only appear when the production coefficients 
are replaced by their complex conjugate for one choice 
of sign in Eq. (IB8I) . That is, if u\ = (b, a\,a 2 ,a 3 ) and 
u 2 = (b',a[,a 2 ,a 3 ), the solutions { 141 , 1 x 2 } and {u±,u 2 } 
should give consistent values for h 3 ± h 2 . This requires 
that either 


h’ 0 + h’ 2 = h* + h* 2 
h’o ~ h 2 = h 0 - h 2 

or 

Hq + h 2 = h 3 + h 2 

h 3 ~ h 2 = h 0 — h 2 

Therefore, either 

3a; -a' 3 = 3a; - a* 3 
b — a l + a 3 = b — a x + a 3 
3 o !^ <3,3 = 3 cii -|- (33 

b f — cii — <33 = b — a 1 — <33 

or 

3 di — (33 = 3 cti — (33 
b' — a\ + a' 3 = b — a\ + a 3 
3ai + a' 3 = 3 a{ + a* 3 
b' — a' ± — a! 3 = b* — a\ — a 3 . 


(BIO) 


(Bll) 


(B12) 


(B13) 


Both Eq. (IB12I) and Eq. (IB13I) require that 


Im b = —2 Im a\. (B14) 

The difference between Eq. (1B12I) and Eq. (IB13I) is a sign 
change for imaginary part of each amplitude. This differ¬ 
ence is equivalent to the trivial ambiguities discussed in 
section llVDI Let us choose the phase convention given by 
Eq. (IB12I) . Finally, invariance of /i, given the conditions 
above, requires that a' 2 = a 2 . 


Using the conditions in Eq. (IB12I) and the constraint 
a 2 = a 2 , the alternate set of solutions can be written in 
terms of the original set as 

Re Vq A = Re Vo.i 

Im V’ A = —^=(3 Im U 2 ,i — Vb Im V 2t2 + Im V 2j3 ) 
Re V 2 l = Re U 2 ,i 

Im V 2 1 = Im V 2 1 + Im V 2 2 + ^ Im V 2 3 
’ 3 o 

Re U 2 ' 2 = Re U 2 ,2 

V5 

Im V 2 2 = —Im U 2 ,2 -— Im U 2 ,3 

Re U 2 ' j3 = Re U 2 , 3 
Im V 23 = Im V 23 . 


(B15) 


Note that the last two lines of Eq. (IB15I) indicate that 
the ambiguous solution for the 2 ++ E3 amplitude is re¬ 
dundant with the original solution. That is, the 2 ++ E3 
amplitude does not exhibit multiple solutions. 

In a practical sense, these results are useful to compare 
the mathematical predictions to what is found experi¬ 
mentally. Essentially, the predicted ambiguous partner 
for a set of fit results in a given bin may be calculated in 
the following way. First, the results must be rotated in 
phase space such that the condition in Eq. (1B14I) is sat¬ 
isfied. Next, the ambiguous partner may be determined 
using Eq. (IB15I) . Finally, this predicted solution must be 
rotated back into the original phase convention. Now, the 
predicted ambiguous partner may be compared with the 
experimentally determined fit results. Studies show that 
the mathematically predicted ambiguities match those 
found experimentally. 


Appendix C: Supplemental Materials 

In addition to the figures presented here, the results 
of the mass independent analysis in each bin of M r o T o 
are included in the supplemental materials [37]. This 
includes the intensities of each amplitude and the three 
phase differences for each bin of o T o. The two ambigu¬ 
ous solutions of the nominal results are separated into 
two text files, while one additional text file contains the 
alternate results in the region where they are not redun¬ 
dant with the nominal results. Note that these results 
contain only statistical errors. 

It is important to reiterate that errors reported in the 
supplemental results (and in the figures in the text) are 
derived from the covariance matrix of the fit parameters. 
That is, they are valid in the Gaussian limit, a limit that 
cannot be guaranteed for all parameters in the analy¬ 
sis. Therefore the use of these results in a subsequent fit 
to parameters of interest cannot be expected to produce 
statistically rigorous values of the parameters. Likewise 
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a x 2 or likelihood-ratio test of a model describing the 
results cannot be rigorously constructed. 

An attempt to quantify the potential systematic bias 
in subsequent analyses was made as follows. First, a sam¬ 
ple of MC with equivalent statistical precision to the data 
was generated using a model consisting of a coherent sum 
of Breit-Wigner resonances in a way that best approxi¬ 
mates the data. A mass independent amplitude analysis 
was performed on this MC sample using the same pro¬ 
cedure that was applied to the actual data reported in 
this analysis. The results of this mass independent anal¬ 
ysis of the MC sample were then fit with a Breit-Wigner 
model, the same model with which they were generated, 
where the couplings of the Breit-Wigner distributions in 
the model were allowed to float as free parameters. While 
most fit parameters exhibited typical Gaussian fluctua¬ 
tions about their known input values, there were some 
non-Gaussian outliers. About one-third of the parame¬ 
ters exhibited deviations from input at or above the three 
sigma level. In comparison with a mass dependent analy¬ 


sis, in which the Breit-Wigner model is directly fit to the 
same mock data, the parameter errors in the model fit to 
the MI results were generally larger, typically within a 
factor of two, but in some cases by up to a factor of ten. 

To probe the scale of the systematic deviations of the 
fitted values from the true input values used to gener¬ 
ate our MC sample, for each amplitude we used the true 
value of the coupling instead of the fitted value and com¬ 
puted (1) the total intensity integrated over all phase 
space and (2) the fit fraction (ratio of individual ampli¬ 
tude intensity to total intensity). We observe the de¬ 
viations in (1) to be at or below the 1% level for all 
amplitudes and deviations in (2) to be at or below 2% 
on an absolute scale for all amplitudes. For small am¬ 
plitudes, this means that relative deviations in intensity 
may occur at a level of 10-90%. This suggests validity 
and precision at a level sufficient for model development; 
however, rigorous values for any model parameters can 
only be reliably obtained by fitting the given model di¬ 
rectly to the data. 
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